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ABSTRACT 



Context. Owing to their computational simplicity, models with elliptical potentials (pseudo-elliptical) are often used in gravitational 
lensing applications, in particular for mass modeling using arcs and for arc statistics. However, these models generally lead to negative 
mass distributions in some regions and to dumbbell- shaped surface density contours for high ellipticities. 

Aims. We revisit the physical limitations of the pseudo-elliptical Navarro-Frenk-White (PNFW) model, focusing on the behavior of 
the mass distribution close to the tangential critical curve, where tangential arcs are expected to be formed. We investigate the shape 
of the mass distribution on this region and the presence of negative convergence. We obtain a mapping from the PNFW to the NFW 
model with elliptical mass distribution (ENFW). We compare the arc cross section for both models, aiming to determine a domain of 
validity for the PNFW model in terms of its mass distribution and for the cross section. 

Method. We defined a figure of merit to i) measure the deviation of the iso-convergence contours of the PNFW model to an elliptical 
shape, ii) assigned an ellipticity £j to these contours, Hi) defined a corresponding iso-convergence contour for the ENFW model. We 
computed the arc cross section using the "infinitesimal circular source approximation". 

Results. We extend previous work by investigating the shape of the mass distribution of the PNFW model for a broad range of 
the potential ellipticity parameter e and characteristic convergence k^. We show that the maximum value of e to avoid dumbbell- 
shaped mass distributions is explicitly dependent on , with higher ellipticities (e = 0.5, i.e., ej ^ 0.65) allowed for small k^. We 
determine a relation between the ellipticity of the mass distribution £ S and e valid for any ellipticity. We also derive the relation of 
characteristic convergences, obtaining a complete mapping from PNFW to ENFW models, and provide fitting formulae for connecting 
the parameters of both models. Using this mapping, the cross sections for both models are compared, setting additional constraints on 
the parameter space of the PNFW model such that it reproduces the ENFW results. We also find that the negative convergence regions 
occur far from the arc formation region and should therefore not be a problem for studies with gravitational arcs. 
Conclusions. We conclude that the PNFW model is well-suited to model an elliptical mass distribution on a larger s—K^ parameter 
space than previously expected. However, if we require the PNFW model to reproduce the arc cross section of the ENFW well, the 
ellipticity is more restricted, particularly for low . The determination of a domain of validity for the PNFW model and the mapping 
to ENFW models could have implications for the use of PNFW models for the inverse modeling of lenses and for fast arc simulations, 
for example. 

Key words. Gravitational lensing: strong - Galaxies: clusters: general - Galaxies: halos - dark matter 
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1. Introduction 

Gravitational arcs are powerful probes of t he mass dis- 
tributi on in galaxies dKoopmans et alJ [2009; Barnabe et al. 
Suvu et all 120121) and galaxy 



2011; 

Mira lda-Escudel Il993ab lHattori et al 

cos mological models dBartelmann et al. 



clust ers dKovnerl lj 
1997) and can be 



to constrain cos mological m odels (IE 
Qguri ~et~aT I200U IGolse et alJ 120021: 1 



_ _ , Bartelmann et al. _ 

Julio et al]l201Q\ The main techniques employed to extract in- 



formation from gravitational arcs have been a rc statistics (i.e. 
counting the number of arcs in lens samples. IWu & Hammerl 
ll993tlGrossman & Sahall 19941: iBartelmann & Weisslll994 and 
inverse modeling (i.e. "deprojecting" the arcs i n individual 
lens systems to determine the lens and the source. [K neib et all 
1993: lKeetonll200ll: IGolse et al.ll2002t IWavth & Webster! 12006: 



Julio et alJ 2007, 2010). The first requires large samples of arcs 



while the second needs detailed information on the lensing sys- 



tems, typically imaging from space and lens and source red- 
shifts. 

These applications hav e triggered arc sea r ches in sur- 
veys covering large areas dGladders et alJ 120031: lEstrada et ail 
2007: ICabanac et alJl2007t iKneib et al. 11201 OtlMore et alJl201 lb 
Belo kurov et al. I 120091: iKubo et al] 1201 Oh. in large spectro - 
cospic surveys dBolton et alJ 20081 : Brownstein et al. [2012I) 
and in sur veys targeting kn own c lusters dLuppino et alJll999t 
Zaritsky & Gonzal ez! 120031 : ISmith et alJ 120051: iHennawi et al 
20081: iMakler et alJ 120167 iRichard et al.M 2010: K ausch et al 
120101 : iFurlanetto et alJ 120121) . Upcoming wide field imaging 
surveys, such as the Dark Energy Survey (Ab bott et al.1 12005: 
lAnnis et ail 120051) and the Large Sy noptic Survey Telescope 
(Stub bs et al.ll2004t Hvezic et alJl2008l) . will lead to the identi- 
fication of larger samples of arcs in thousands of galaxies and 
galaxy clusters, well suited for arc statistics. Moreover, deep ob- 
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servations from space combined with massive spectroscopy have 
been obtained for a limited number of clusters and were used 
for detailed mass modeling (see, e.g ., iBroadhurst et alJ feOQS; 
iRichard et al.ll2010H.Tullo et al.ll2010h . 

The simplest models that can account for some observed 
properties of arcs (multiplicity, relative positions, morphology) 
are built from axial models by introduci ng an ellipticity either on 
the mass distribution (elliptical models. l Schramm 1990; Ba rkanal 
Il998t iKeetonl 120011: lOeuri et all [20031) or on the lensing po- 
tential (pseudo-ellitpical models , [Blandford & Kocha nek 1987; 
Kassiola & Kovner 1993; Kneib 120011) . Most parametric an aly- 
ses of arcs, both for the inverse modeling (|jullo et al.ll2007l) and 
for arc statistics dOgurill2002l ; lOguri et afl l2003). involve one or 
more elliptical/pseudo-elliptical models (adding, in some cases, 
external shear and substructures). 

Elliptical models, whose surface density is constant over el- 
lipses, are more realistic than pseudo-elliptical ones. For ex- 
ample, elliptical models are motivated by the results of N- 
body simulations, which show that dark mat ter halos are tri- 
axial dJing & Sutol 120021; iMaccio et ail 120071) . such that their 
overall mass distribution can be modeled at first order by ellip- 
soids, whose surface density contours are elliptical. In contrast, 
the surface density of pseudo-elliptical models generally has a 
patholo gical behavior, exhibiting re gions where it takes negative 
values ( Blandf ord & Kochanek! 19871) and prese nting a "dumb- 
bell" (or "peanut") shape for high ellipticities (lKovnerl ll989; 
ISchneider et alJl992t lKassiola & Kovnerll99"3l which does not 
represent the mass distribution of most physical systems. 

On the other hand, pseudo-elliptical models provide sim- 
ple analytic solutions for some lensing quantities, allowing for 
fast numerical methods to be implemented, whereas elliptical 
models require the com putation of integrals, which a re more de- 
manding numerically (Schra mm! 1990tlKeetonl200 lb . Therefore, 
studies that require numerous evaluations of the lensing quanti- 
ties often employ pseudo-elliptical models. For example , sev- 
eral s t udies using arc simulations ha ve used these models dOguril 
l2002t iMeneghetti et alJl2003L 2007 ). Popular codes for lens in- 
versio n dGolse et alj|2002t IWavth & Webster! l2006t Uullo et alj 
120071) are implemented using this type of model, too. 

It is therefore relevant to determine a "domain of validity" 
for pseudo-elliptical models such that the negative convergence 
appears far from the arc-forming region and the shape of their 
mass distribution is closer to elliptical. The determination of 
such a validity region could be useful, for example, to evaluate if 
a set of model parameters derived from the inversion of a system 
with arcs using the PNFW model is physically ac ceptable. Even 
for methods that u se the ENFW for lens inversion dKeetonl200lt 
I Suvu et alJl20H . the PNFW could be useful, within its domain 
of validity, for a faster coarser probing of the parameter space 
that would subsequently need to be refined with the elliptical 
model. 

Once a pseudo-elliptical model is found to be acceptable, it 
is nevertheless necessary to provide a correspondence to an el- 
liptical model. This would be necessary, for example, to compare 
results derived from the mass modeling using observed arcs with 
theoretical predictions. We therefore need to establish a mapping 
among the parameters of the two models. This mapping could 
also be used to replace an elliptic model by its corresponding 
pseudo-elliptic in lensing simulations. 

In this paper we focus on the widely used Pseudo-Elliptical 
Navarro-Frenk- White (hereafter PNFW) model and investigate 
its domain of validity as well as the mapping to the correspond- 
ing elliptical model (hereafter ENFW). To test the equivalence 
among the two models we compute the cross section for arc 



formation, from which we obtain additional constraints on the 
model parameters. 

The outline of this paper is as follows: In Sect. [2] we briefly 
review the PNFW lens model, introducing the conventions and 
parameter ranges to be used throughout this work. In Sect. [3] we 
discuss the region of arc formation and derive physical limits of 
the PNFW mass distribution in this region. In Sect. [4] we con- 
sider two ways for assigning an ellipticity to the PNFW surface 
density. In Sect. [5] we obtain a mapping from the PNFW to the 
ENFW models. In Sect. [6] we compare the arc cross sections 
of the two models using the mapping relations. In Sect. [7] we 
present the summary and concluding remarks. In Appendix lAl 
we present useful relations to derive some lensing functions for 
pseudo-elliptical models. In Appendix|B]we provide fitting for- 
mulae for the limits on the mass distribution of the PNFW model 
and for the mapping to the ENFW. 

2. The pseudo-elliptical NFW lens model 

2.1. Circular NFW lens model 

From N-body simulations Navarro, Frenk & White (NFW) 
found that the radial (i.e. angle-averaged) density profile of 
dark matter haloes approx imately follows the universal function 



dark matter naloes approxin 
dNavarroet al.ll 1996L 1 1997h 



P(r) = 



(r/r s )(l+r/r s ) 



2' 



(1) 



where r s is the scale radius and p s is the characteristic density 
of the halo. This profile has been widely used to represent the 
dark matter mass distribution galaxy to cluster scales and will be 
employed throughout this work. 

From the density profile dTJ, defining the dimensionless 
Cartesian coordinates in the lens planeQ x = £/r s , the conver- 
gence, deflection a ngle, shear, and lensing potential are given by 
(Bartelman3l996|): 



k(x) — 2k s F(x), 

t ^ a 
a(x) — 4k s , 



y(x) = 2 k, 



2g(x) 



-F(x)\, 



ifi(x) — 2k s Ii(x), 



(2) 
(3) 

(4) 
(5) 



where the functions F(x), g(x) and h(x) are defined in 
iGolse & Kneibl d2002l hereafter GK02) and the characteristic 
convergence, k s , is given by 



p s r s 

K s = =— , 
•^cnt 

with the critical surface mass density 
- 2 D ( 



■^-nrtt — 



>OL 



(6) 



(7) 



where Dol, Dls and Dos are the angular-diameter distances 
between the observer and lens (at redshift zl), lens and 
source (at redshift zs), and observer and source , respectively 
(ISchneider etal .1119921: fMollerac h & Rouletll2002h . 



1 In Eq. {D r = V?" 

+ z 2 , where z is the coordinate on the lens- 
observer direction and f = |f | is the coordinate perpendicular to the 
line-of-sight (lens plane). 
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Notice that with the choice of dimensionless coordinates x, 
the lensing functions are independent of r s . Naturally, dimen- 
sional quantities have to be scaled accordingly to recover their 
physical units. For example, angular quantities have to be multi- 
plied by r JDoL to be in radians. 

2.2. Pseudo-elliptical models 

To construct "elliptical models" from a given radial density pro- 
file, the radial coordinate x is replaced by 



9 9 

q + «2 xi, 



I 2 

such that the ellipticity is given by 



(8) 



(9) 



where we assume that 02 > fli> such that the major axis of the 
ellipse is along x\. 

This substitution can be performed on the projected mass 
distri bution, k(x) — » k(x k ), leading to elliptical mass distribu- 
tions dBourassa et aljfl973t iBourassa & Kantowskilll975l: iBravl 
11984 . An alternative is to construct pseudo-elliptical models, by 
making the substitution in the lensing potential, <p(x) — * <p(x e ), 
which leads to simple analytic expressions for the lensing func- 
tions (GK02, iKassiola & Kovnerll 19931: IOgurill2002t iJullo et al.1 
120071) . In particular the convergence and the components of the 
shear are given by (see AppendixlAl 



K e (x) = JIk(x e ) -Sy(x E )cos2(f> B , 
yi s (x) = SK(x E )-^ly(x s )cos2(p s , 

y lE (x) = -^m 2 -S 2 y(x s ) sin20 e , 
yl(x) = M 2 y 2 (x e )-2M!BK(x e )y(x e )cos: 
+S 2 [K 2 (x E )-sm 2 2<p e y 2 (x B )], 



(10) 
(11) 
(12) 

(13) 



-(«! - 02), and k{x e ) and y(x s ) are 



where 3i — \(a\ + a.2), !L - - 
the convergence and shear of a circular model evaluated at x e , 
respectively. 

The PNFW model is obtained by introducing the ellipticity 
on the NFW lens potential (Eq. [5}, where we denote the charac- 
teristic convergence (Eq.[6]l by k%. 

In this paper w e adopt the convention 
(Blandfor d"&Kochanekl[ l987. GK02) 



fll = 1 — S, fl2 — 1 + B, 



(14) 



such that the ellipticity of the lensing potential (Eq. |9j is given 
by 



1 



1 +e 



(15) 



and the ellipticity parameter e is defined in the range < s < I. 

Using this convention, Eqs. (TT~Qb and dT3l > yield Eqs. (17) and 
(19) of GKO20 Another common choice of parameterization is 
a\ = 1/(1 - s), «2 = 1 — e. In this case the ellipticity of the 
potential is simply given by b v — s and Eqs. (1 1 OH 1 31 > yield the 



2 After accounting for the known typo, in Eq. (19) of 
GK02 (cos 2 20 £ — > sin 2 2(p s ). We tha nk the a u thors for 
pointing this out on the lenstool code iJullo et aT] l2007t) at 
http://www.oamp.fr/cosmology/lenstool/ 



same expressions for the lensing functions as, for example, in 
iLima etaU(l201Ctfl 

To obtain the maximum value of to be used in this pa- 
per we considered the most extreme lensing systems expected 
in nature, wiftfl M 2 oo = 4 x 10 I5 /i _I M o , zl = 16, Zs =7 
(Oguri & Blandford 2009). The values of r s were obtained from 
the distribution of the concentration parameter c = r^nn/r?. 
p(c\M , z = 0), derived from N-bod y simulations by Neto et all 



(2007) with redshift scaling given in Maccio et al. (2008). From 
M200, c, zl, and zs we obtained the characteristic convergence 
by applying the rel ations for the sp herical NFW model (Eqs.Q] 
|5] and |7] see, e.g., ICaminha et al. I l20l2h . where we chose the 
ACDM matter and cosmological constant density parameters as 
Q m = 0.3 and Qa = 0.7, respectively. We generated many real- 
izations of the C-M200 relation, converted them to characteristic 
convergence, and took the 95% upper limit of the derived dis- 
tribution as its maximum value, which leads to k £ =a 1.5. This 
sets an inclusive upper limit for k v s , at least within the ACDM 
framework. 



3. Physical limits of the PNFW mass distribution 

We focused on the mass distribution in the vicinity of the re- 
gion of tangential arc0 formation. The arcs are usually de- 
fined as an image with length-to- width ratio L/W above a 
given threshold R±. For infinitesimal circular sources this ra- 
tio can be determined from the radial and tangential eigenval- 
ues of the J acobian matrix of the lens mapping, A r and A t , 
respectively dWu & Hammerlll993t iBartelmann & Weissl ll994: 
Ham ana & Futamasdll997 ) : 



— = 1^1 
W 



(16) 



where :- A r /A,, with A r = 1 — K + y and A t = 1 — K — y. 

Fixing a value for R^ determines a region limited by the 
curves R^ = +R t h (constant distortion curves), where gravita- 
tional arcs are expected to be formed. Although condition (fT~6b 
does not hold for images of sources crossing the t angential caus- 
tic (merger arcs, Ro zo et al. 2008; F erreir a 2010), nor for large 
or noncircular sources, the curves defined above still provide a 
typical scale for the region of arc formation. A common choice 
for the threshold is R t h = 10, which we adopted in this work (un- 
less explicitly stated otherwise). The tangential critical curve is 
given by the condition R^ = 00. In Fig.[TJthese curves are shown 
for a few combinations of k% and e. 

Once a domain for arc formation has been established, we 
need to associate to this region iso-convergence contours, de- 
noted by k e contours, which define the shape of the mass dis- 
tribution. We chose to match the k e contours to the R^ = ±R± 
curves at the major axis (X2 = 0), since most a rcs are ex- 
pected to be formed close to this region (see, e.g., Dalai et al 



2004 lComerford et al.l2006HMeneghetti et al.l2007trMore et al 



2011). This definition is illustrated in Fig. Q] where distortion 
curves and their associated iso-convergence contours are shown. 



3 After rotating the lens by n/2 to follow their convention. 

4 Where M 2 oo is defined as the mass contained within a radius r 2 oo 
enclosing a region with mean density 200 times the critical density of 
the Universe at zl and h is the Hubble parameter in units of 100 Mpc 
krrr's. 

5 Radial arcs are more difficult to observe because they are hidden 
by the light of the lens, since they are forme d in the central region 
of the lenses and are usually fainter images ( Miralda-Escude 19911 : 
IBarte lmann 200j). 
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Fig. 1. Critical curve (R^ = oo), curves of constant distortion (R \ = ±/?th) and k e contours associated to each curve (solid lines) 
for R± = 10 and = 0.1 and s = 0.45 (left panel), «f = 0.8 and s = 0.35 (middle panel) and - 1.5 and s = 0.3 (right panel). 
The axes are in units of r s . 
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Fig. 2. Negative convergence of the PNFW model. Solid lines show the intersection of k e — contours with X2 axis. Other lines 
correspond to the intersection of the iso-convergence contours associated to R^ - R& curves with the X2 axis. Left panel: k s = 0.1. 
Right panel: k s = 1.5. 



We will use this choice throughout this paper and refer to the k e 
contour by the value of associated to it. 

A fundamental problem of pseudo-elliptical models is 
the presence of region s with negative mass distribution 
dKassiola & Kovnerlll993h . For the PNFW model parameterized 
as in Eq. (TBI , negative k e contours form lobes oriented along the 
X2 axis. These regions occur far from the lens center and for any 
s > 0. The k e — contours are independent of k% and the nearest 
point to the region of tangential arc formation is located on the 
X2 axis. For example, in Fig. Q] for s = 0.3, 0.35, and 0.45, neg- 
ative values of k e arise at xi = +10.5,7.3 and 4.3, respectively, 
well outside the range of these plots. 

We have verified that the k e — contours do not intersect 
the iso-convergence contours associated to the \R,\\ = R^ curves 
in the whole e and range, as can be seen on Fig. [2] As ex- 
pected, the negative convergence lobes approach the R^ = const. 
curves as the ellipticity increases, but never intersects these re- 
gions, even down to R t ^ — 1.25. The lower the characteristic 



convergence, the farther is the negative convergence region from 
the arc formation region. Therefore, the formation of tangential 
arcs occurs far from the k e < regions, and can be ignored in 
the context of this paper. 

Now we turn to the shape of the k e contours. Given the cho- 
sen orientation of the major axis, for contours that are close to 
elliptical, the maximum value of X2 (x™ ax ) is located at X\ = 0. 
On the other hand, for dumbbell-shaped contours x™ x is located 
at xi + 0. This simple property can be used to determine the 
maximum value of s (e max ) such that for s > e max a dumbbell- 
shape emerges. This is shown schematically in Fig. [3] For a given 
Kg, we started with a low value of s and computed x™ ax for the k e 
contours associated to each R\ curve. As the ellipticity parame- 
ter is increased, e max is attained when the point corresponding to 
x™ ax starts to be located at x\ + 0. Repeating this procedure for 
any value of k% and R\ leads to the function e max (^i , Ra)- 
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We found that this function is weakly dependent on Rj. 
For example, the absolute difference \s max {i<i! , R A = co) — 
e max (/Cs,/? / i = ±Rth)\ is at most 0.03 for R t h = 4, corresponding 
to a fractional difference of about 10%, in the whole k% range. 
This maximum difference decreases to 0.01 for R± = 10. Thus, 
it is sufficient to define £ max (Kj) '■= e max (K^,Rji = oo) in the arc 
formation region. This function is shown in Fig. [4] Clearly e max 
depends on k v s , decreasing for high characteristic convergences. 
For low values of k v s , e max converges to 0.5. In Appendix lB.il we 
provide a best-fitting function to E mas («-j), which could be use- 
ful, for example, to check if a given solution from inverse mod- 
eling has a physically meaningful mass distribution. The limits 
on the shape of the iso-convergence contours will be revisited in 
the next section by measuring the deviation with respect to an 
elliptical shape. 



Rj. = R 

R,=-R. 



"th 



R, = oo 

R, = -R, 
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0.6 
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D 2 = 0.45 x 10" 3 
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0.25 0.5 



0.75 
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Fig. 3. Sketch of the method for determining the maximum 
value of s to avoid dumbell-shaped k s contours associated to the 
R,i = const, curves. Left panel: Shape of the iso-convergence 
contours for e < e max . Right panel: Shape of the iso-convergence 
contours for s > s lmx . 



4. Ellipticity of the PNFW mass distribution 

We adopted two procedures to associate an ellipticity to each 
k b contour and to measure its deviation from an elliptical shape. 
In the first, we followed GK02 and define the semi-major axis 
gegk and semi-minor axis Z?gk as the intersections of the k e con- 
tour with the x\ and X2 axes, respectively (see Fig. |5), such that 
the ellipticity is 



OK 



1 



bGK 
flGK 



(17) 



The second procedure is to fit the k s contour by an ellipse. 
For this sake we introduce a figure-of-merit that represents the 
mean weighted squared fractional radial difference between the 
contour and the ellipse 

2 ._ Zli wtirm - rsM? 



D l := 



where N is the number of points on the k s contour, fa is their 
polar angle, w,- = fa - fa-\ is a weight accounting for a possible 
non-uniform distribution of fa, r((pi) is the radial coordinate of 
the k e contour and r%(fa) is the radial coordinate of the ellipse, 
given by 



rs = 





1 sin fa 


2 




V b J 





-1/2 



(19) 



Fig. 4. Maximum value of the ellipticity to avoid dumbbell- 
shaped mass distributions in the region close to the R A = oo 
contour as a function of k%. The solid line corresponds to e max 
obtained from the procedure in Sect. [3] The dashed line cor- 
responds to the values of e obtained from upper limits of the 
figure-of-merit D 2 representing the fractional deviation of the k b 
contours with respect to an ellipse (Sect.|4j>. 



where a and b are the semi major and minor axes, respectively. 
The form of D 2 in Eq. ( fT8b was chosen to be scale-invariant and 
independent of the discretizatiorQ. Owing to the symmetry of the 
K e contour, it is sufficient to compute D 2 in the first quadrant of 
the lens plane. 

The best-fitting ellipse is found by mini mizing D 2 , for which 
we used the MINUIT code (Jame s 1998). The resulting best- 
fitting values from this elliptical fit (EF) method, aef and b^e, 
yield the ellipticity 



1 



<zef' 



(20) 



When the iso-convergence contours are close to elliptical, 
we expect the results from both methods to be very similar. 
However, differences could emerge for high values of e, espe- 
cially in the dumbbell-shape regime. Another difference that will 
be discussed in Sect.|5]arises when assigning a k e contour to the 
derived ellipse. 

We may associate an ellipticity to the iso-convergence 
contours related to each distortion curve. However, the func- 
(18) tions £ £ ' (s,Ks,Ra) are weakly dependent on R A close to 
the arc formation region. For example, the absolute differences 
\e^ K ' EF (e,KtRA = c») - e£ K ' EF (e, /rf , ±R tb )\ are at most 0.01 for 
R t h = 10 and 0.02 for R± = 4. Therefore, we chose the ellipticity 
associated to the convergence at the critical curve (R^ = oo) as 
the ellipticity in the arc formation region. In Fig.|6]the resulting 
values for e^ K (dashed lines) and s EF (solid lines) are shown as 
a function of e, for «f = 0.1, 0.5, 1.0, and 1.5. 



6 A convergence to within about l c , 
parameter ranges considered here. 



is achieved for N = 100 for the 
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Fig. 6. Ellipticity of the PNFW mass distribution from the elliptical fit (ef F , solid lines) and GK {s^ K , dashed lines) methods, as 
a function of the parameter s for four values of k%, calculated at the tangential critical curve (R \ = oo) at = 0. The dotted line in 
the upper left panel corresponds to the ellipticity of the lensing potential Sy{s), Eq. ( IT5b . The dot-dot-dashed line shows the ey. — s 
line to guide the eye. The vertical dot-dashed lines show the values of s max (K%). 



As expected, the two ellipticity measures agree very well for 
low values of s. For example, for s = 0.25, |e FF - * s at most 
0.03 on the whole /c, range. As can be seen in Fig. [6] although the 
behavior of the two functions is qualitatively very similar in the 
whole ellipticity range, noticeable differences appear precisely 
for the values of e (shown in Fig.|4j> close to where the dumbbell- 
shape arises. 

The minimum value of the figure-of-merit ( fT8b £^, in can be 
used as a goodness of fit and hence as an estimator of the depar- 



ture from the elliptical shape. By setting a threshold on T> ^ in 
a maximum value of s can be determined, for each k v s , such 
that D 2 does not exceed this threshold, ensuring a small devia- 

min ' to 

tion from the elliptical shape. In particular, setting this threshold 
at 4.5 x 10~ 4 avoids the dumbbell-shaped mass distribution, as 
shown in Fig|4] Thus, a small deviation from the elliptical shape 
is deeply connected to the avoidance of dumbbell shapes. Both 
conditions impose similar restrictions on the ellipticity and im- 
ply that the PNFW could be used to model the mass distribution 
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EF 



k contour 

e 

GK method 
EF method 



GK 




a GK \ \ a EF 



Fig. 5. Illustration of the methods used to associate the PNFW 
mass distribution to elliptical contours. In the GK method, aqk 
and £>gk correspond to the semi-major and semi-minor axes of 
the ellipse (dashed line) passing through the intersection of the 
iso-convergence contour (solid line) with the X\ and X2 axes. In 
the EF method, oef and Z?ef correspond to the semi-major and 
semi-minor axes of the best-fitting ellipse (dash-dotted line) ob- 
tained by minimizing Eq. ( fl~8b . 



in the region of arc formation for s below the limits given in 
Fig.l4land Eq. (HTHi. 

The validity of the PNFW to represent elliptical mass distri- 
butions was also investigated in GK02. They considered a sin- 
gle value of the characteristic convergence and investigated the 
shape of the mass distribution as a function of the distance of the 
k e contour to the lens center. They provide a fitting function for 

4 



a GK + ^gk anc ^ e ' va lid f° r s < 0.25. 



sy. as a function of xqk 
Their lensing potential depth is expressed in terms of a charac 
teristic velocity v c , connected to the NFW profile parameters by 
(iGolse et alJl2002h 



(21) 



Using the values in GK02 (v c = 2000 km/s, r s = 150 kpc, zl — 
0.3, and zs = 1), assuming Q m = 0.3, Q.\ = 0.7, and Ho = 
65 km/s/Mpc, and using Eqs. © and ^ yields k% - 0.88. Our 
results for e^ K (e, - 0.88) reproduce their fit to within 2% 
down to its limit of validity. 

On the other hand, by exploring a large interval of *f , we 
find that higher values of s may be allowed, at least in the region 
of arc formation. In particular, values as high as e 0.5 (corre- 
sponding to bj_ =t 0.65, see Fig.|6]l are permitted for low k v s . It is 
therefore useful to provide a fit that is valid beyond s = 0.25 and 
includes the dependence of the ellipticity with *f . Such a fitting 
function for s^(s, K t) is shown in Appendix IB. 21 which is valid 
in the whole range of s and the range of *f considered in this 
paper. 



5. Mapping among the PNFW and ENFW models 

The ENFW model is constructed by replacing k (x) in Eq . [2] by 
k(x^) (see Sect. 12.2b . 
2012) 



where x ei . is chosen as (Caminha et al. 



(1 -e£)x\ + 



1-es 



(22) 



such that ex is the ellitpicity of the mass distribution. In this case 
we denote the NFW characteristic convergence by a^. 

We constructed a mapping among the PNFW and ENFW 
models such that their mass distribution is similar on the arc for- 
mation region. In other words, for each pair (k%, s) we associ- 
ated a corresponding pair of the ENFW model parameters (a^, 
sz), where a^ is the characteristic convergence of the associated 
ENFW model. The ellipticity of the ENFW mass distribution 
is simply given by the ellipticity associated to the k e contours, 
as discussed in Sect. [4] with approximate expressions given in 
Appendix IB. 21 

The determination of *f is numerically more complex and 
involves an ambiguity on how to perform the matching among 
models. We have chosen to perform the matching at oef, i. e., at 
the intersection of the best-fitting ellipse to the k s contour asso- 
ciated to a given R A = const, curve with the x\ axis. We have 
considered two possibilities: i) matching the value of the ENFW 
convergence at «ef with the PNFW convergence associated to 
the = const, curve, if) matching the position of = const. 
curve of the ENFW model, i.e. such that it intersects the x\ axis 
at aef- For any pair (**, s) we fixed the ellipticity of the ENFW 
model as sz(s, a^) and obtained *f following the two procedures 
above. 

The two possibilities described above yield very similar val- 
ues for a^ down to high ellipticites. For example, taking as ref- 
erence the k e associated to the critical curve (i.e. to = °°), 
the maximum difference of among the two procedures for 

= 1.5 and e = 0.8 is 3% and this difference decreases sub- 
stantially for lower ellipticities and k^. On the other hand, if in- 
stead of using «ef as reference position for the association of 
we use qgk, the results are still very similar for s < e max (^), but 
may differ substantially for higher ellipticities. 

As described above, convergences can be matched for k s 
contours associated to any R^ curve. We verified that the de- 
rived characteristic convergence is almost constant in the re- 
gion of tangential arc formation. Indeed, the absolute differences 
\k%(s,ic*,R a = oo)- ^(£,4",^ = +R th )\ are at most 0.005 (0.01), 
for < 0.05 and 0.01 (0.03), for *f ~ 1.5, with R th = 10 
(R t h = 4), and e < 0.6. Therefore, as in the case, we chose the 
characteristic convergence in the region of arc formation as the 
value of *f {s, k%) calculated at the intersection of R A = oo with 
the x\ axis. As will be discussed in Sect. [6] we chose method (if) 
to make the association between the two models. 

In Fig. |7] we show as a function of a* for some values 
of e. As expected, for e — we have *f = *f , but this equality 
does not hold for non-zero ellipticities. In particular, is always 
larger than its corresponding <cj, and the difference between them 
increases with s. In Appendix IB. 3 1 we provide a fitting function 

for Aj (fi, Aj). 



6. Comparison of the arc cross section 

The efficiency of a lens to produce arcs is quantified by the arc 
cross section <x Rlh , which is defined as the area in the source plane 
that generates images with L/W > R t h, weighted by the multi- 
plicity of the images (i.e. multiply imaged regions are counted 
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Fig. 7. Relation among the characteristic convergences of the 
PNFW and ENFW models for some values of the parameter e 
calculated at the tangential critical curve (Rj = oo) at X2 = 0. 
The solid line shows the /cf = k? line to guide the eye. 



multiple times, see e.g.. iMeneghetti et alj |2003). The computa- 
tion of the arc cross section usually requires e xtensive arc sim- 
ulations, which are computationally expensi ve (Miralda-Escude 
1993 btlBartelmann & Weissll994tlMeneghetti et alj200 11.120031; 
Qguri et al]l2003h . However, in the infinitesimal circular source 
approximation, Eq. dT6t . the cross section can be obtained di- 
rectly from the local mapping from lens to source plane. In this 
case &-R th is easily computed in the lens plane (in dimension- 
less c oordinates) by (see, e.g. jFedeli et alj|20"pl^ ICaminha ( 
l20T2h 



L etalJ 



HxT l d 2 x, 



(23) 



where ju = (A r x A,)~ l is the magnification and the integral is 
performed over a region with local distortion above the given 
threshold (i.e., the arc formation region introduced in Sect.[3]and 
used in the preceding sections). Since the mapping between the 
PNFW and ENFW models was constructed in that region, the 
cross section is well suited to check if this mapping, obtained 
from matching the mass distribution, also holds for other lensing 
quantities, in this case the magnification. 

If the PNFW can be used to replace the ENFW in some appli- 
cations (for example, arc statistics), we would expect the predic- 
tions for <x« lh to be similar for both models, at least in the region 
where the PNFW provides an adequate description for the mass 
distribution. On the other hand, if the predictions do not match 
in this region, they can be used to set additional constraints on 
the PNFW model. 

To compare the cross section for both models we defined 
a regular grid in the (<c? , s) parameter space and mapped each 
point to (/cf , s£) using the procedures described in section [5] 
The cross sections <tpnfw and o"enfw were computed for each 
set of parameters from Eq. ( l23l using the expressions for A r (x) 



and A,(x) of the correspondin g model. The computa tion for the 
ENFW model was taken from Caminha et a i1 (l20T2h . The result 
for both cross sections is shown in the left panel of Fig. [8] where 
contours of constant &R lh are displayed. Visually, the arc cross 
sections of both models are similar in the region s < e max - 

To quantify the difference between ctenfw and <5"pnfw, we 
computed the relative difference 



Ao- 



O" ENFW _ CT PNFW 



C ENFW 



(24) 



We have computed this fractional difference using the var- 
ious possibilities for associating ENFW model parameters to 
PNFW parameters discussed in section [5] Although those def- 
initions generally yield similar values for 1^, the value of A&I& 
can vary substantially, since the cross section is very sensitive to 
k s . Comparing the matches at «ef and «gk, the latter gives a frac- 
tional difference at least 50% higher than the former and this dif- 
ference increases substantially with a - ?. The difference between 
the methods (/) and (ii) is much smaller and is not very sensitive 
to /<:? , but in general procedure (ii) leads to smaller differences 
in the cross section than method (/). We chose therefore to de- 
fine the matching among models using method (if) at «ef and 
compared the cross sections using this choice. 

In the right panel of Fig. |8]the contours of constant relative 
difference are shown. Our results show that for s < e max , A&I& 
can be as high as 30% for low values of a:? . Thus, even in the re- 
gion where the ENFW and PNFW mass distributions are similar, 
there can be substantial deviations on the cross section. 

We may combine the constraints from the shape of the mass 
distribution and by assuming a maximum fractional deviation 
for the cross section. For example, for a maximum deviation of 
10%, the joint constraint leads to a region limited approximately 
by the lines 



0.25 - 1.30k? 



s = { 0.08 + 0.42k?, 0.1 < /c? < 0.65, 



<0.1 



(25) 



0.41 - 0.09k? 



> 0.65. 



Within this region the PNFW model can reproduce both the local 
mapping and the mass distribution of the ENFW model. 



7. Summary and concluding remarks 

Motivated by its potential applications for gravitational arcs, we 
revisited the PNFW model, seeking to determine domains of va- 
lidity in terms of the mass distribution and the arc cross section. 
We have shown that the lensing functions of pseudo-elliptical 
models have simple analytic expressions (Eqs. [T0HT31 for any 
choice of the parameterization of the ellipticity. 

We analyzed the PNFW mass distribution, in the "arc for- 
mation region" limited by the constant distortion curves (R^ = 
±R t h) by associating K e contours to these curves. We verified that 
the convergence never takes negative values in this region. Since 
the results obtained in this work are weakly dependent on R^, 
for Rfh > 4, we chose the critical curve (R^ — > 00) to derive the 
final relations summarized below. 

We determined the maximum value of the potential elliptic- 
ity parameter, e max ('<'?), such as to avoid dumbbell-shaped mass 
distributions. The results (Fig. |4j enlarge the domain of applica- 
bility of the PNFW model to describe the mass distribution, at 
least in the arc formation region, beyond the commonly adopted 
value of e ^ 0.25, allowing values as high as e = 0.5 for low 
values of /c?. 
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Fig. 8. Comparison between arc cross sections. Left panel: Contours of constant arc cross section in terms of the PNFW parameters. 
Solid lines correspond to ixenfw and dashed lines correspond to oenfw- The contours of constant arc cross section from left to right 
are 1 x 10 6 , 1 x 1(T\ 1 x 10~ 4 , 10~ 3 , 0.01, 0.02, 0.04, 0.08, 0.16 and 0.24. Right panel: Contours of A<x/cr (Eq.|24]l. In both plots, 
the dash-dotted line shows SmaxOfj)- Calculations were made for/? t h = 10. 



We introduced the figure-of-merit D 2 (Eq. [T8l to quantify 
the deviation of the k e contours from the elliptical shape. Setting 
a maximum value for D 2 at 4.5 x 10~ 4 also avoids dumbbell- 
shaped mass distribution (Fig @}, showing that the contours de- 
viate from the elliptical shape on the verge of the emergence of 
the dumbbell shape. 

The function D 2 can be used to assign a best-fitting ellipse to 
the k e contour and hence to obtain an ellipticity ej_ for the mass 
distribution (EF method). We showed that the ellipticites ej; ob- 
tained from this method are almost identical to the method in 
GK02, especially for low ellipticities. However, the EF is bet- 
ter suited to assign an iso-convergence contour to the k con- 
tour. Furthermore, using the EF allows one to match between 
the PNFW and the ENFW models in a way that minimizes the 
difference between cross sections. 

We provided fitting functions for e^ F (e, k%) (Eqs. IB.2J4B.3l l. 
extending the results of GK02, in the arc formation region, for 
any value of s and including the dependence on k*. Going to 
higher ellipticities is relevant, given that values as high as s 0.5 
(corresponding to sz 0.65, see Fig. |6]l are allowed in the arc 
formation region. 

From N-body simulations, it is found that the probability dis- 
tribution of the projected ellipticity peaks at about g£ _ 0.5, with 
only a small fraction of the halos having sj_ > 0.6 (Qgur ieTaT] 
2003). Converting e ma x(^) to we found that values of > 
0.5 are allowed in the whole investigated range of k% . Therefore, 
the PNFW would provide a good description of the ENFW mass 
distribution for most expected values of s%. 

By associating the iso-convergence contours of the PNFW 
to the ENFW close to the tangential critical curve, we obtained 
a relation among characteristic convergences k^(s, k%). Fitting 
functions for this relation are provided in Appendix IB . 3 1 (for the 
case (ii) discussed in Sect.|5j- Combined with the relation among 
ellipticities, this function completes the mapping of the parame- 



ters (e, /fj) of the PNFW model to the parameters (ex, ) of the 
ENFW model. 

To test the mapping in a practical application we compared 
the predictions for a quantity that is useful in arc statistics. We 
computed the arc cross sections for both the PNFW and ENFW 
models and compared their predictions by matching the model 
parameters using this mapping. We did not find a direct con- 
nection between Acr/a - and e max (which has a similar shape as 
a function of k% as the contours of constant D 2 ), although the 
cross section was computed in the region where the mapping is 
obtained. In other words, the limits derived from the shape of the 
mass distribution do not match those from cross section. 

We may use this result to set additional constraints on the 
parameters of the PNFW model, by requiring an agreement with 
the ENFW for & in addition to the condition e < e max (k%). This 
would ensure that the mass distribution as well as the local map- 
ping (represented by the magnification /j) are well reproduced 
by the PNFW. Approximate limits of this combined restriction, 
imposing an agreement of about 10% for the cross sections, are 
given in Eq. $25[ . 

This new restriction should now be tested in other applica- 
tions, especially with simulations using finite sources, to check 
if the PNFW and ENFW can be mapped to reproduce the same 
physical results. This would validate the use of pseudo-elliptical 
models for simulations and the inverse problem, providing the 
relation to the associated elliptical model. 
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Appendix A: Lensing functions for pseudo-elliptical 
models 

To derive the lensing functions for pseudo-elliptic models, it is 
useful to introduce the following coordinate transformation: 



X] 



COS ( 



X2 



V«2 



sine 



(A.l) 



where x E = -^ai x 2 x + a 2 x\ and <p £ - arctan(jt2/*i y[a 2 ja\). The 



Jacobian matrix of this transformation is 
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stn( 
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(A.2) 



Since the gradient operator transforms as V x 
^T~ l (x,x s )V x , we have 



d Xl - ^cos(p s d Xe sin0 e 3 0E , 

dx 2 — V^2 sin (p £ d XE H cos <p £ d<p E ■ 



(A3) 
(A.4) 



Using the expressions above, the deflection angle (a s (x) 
V x ip(x s )) for an elliptical potential reads (GK02) 



a e (x) 



u\{x) 
a 2 {x) 



a(x B ) y/aicos (f) e 
a(x e ) -\Jci2 sin (p e 



(A.5) 



where a(x E ) is the deflection angle of a circular model evaluated 
at x — x £ . 

Taking the partial derivatives of each component of the an- 
gle deflection and using Eqs. dA.31 — lA~4l > it is straightforward to 
obtain 



d x .ai(x) = ai 



d X2 a 2 (x) = a 2 



da(x e ) 2 a(x e ) . 2 
cos cf> e H sin <p t 

C$X £ 

da(x e ) 



. 2 , a ( x e) 2 

sin d> F h cos ( 



d X2 ai{x) 



da(x e ) a(x s ) 



sin20 £ 



= d x .a 2 (x). 



(A.6) 
(A.7) 

(A.8) 



Using the relations for the lensing functions for circular poten- 
tials 



<x) = \ 



a(x) da(x) 
x dx 



7(x) = \ 



a(x) da(x) 
x dx 



(A.9) 
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it is possible to express Eqs. ( IA.6I ) 
and y(x E ), i.e. 



A. 8b as a function of k(x b ) 



d Xl ai(x) = a\ [k(x e ) - y(x £ ) cos 20 E ] , 
d X2 0!2(x) = a 2 [k(x s ) + y(x s ) cos 2<f> s ] , 
d Xl a\{x) = - y]a x a 2 y{x s ) sin20 e = d Xl a 2 (x). 



(A.10) 
(A.ll) 
(A. 12) 



Applying the usual definitions for the convergence and shear in 
terms of the deflection angle, we obtain Eqs. (1 1 OH 1 31 >- 



Appendix B: Fitting formulae for e max and for the 
PNFW— ENFW mapping 

B. 1 . Fitting formula for e max . 

Applying the procedure outlined in Sect. [3] we obtained the max- 
imum value e max to avoid the dumbbell-shaped mass distribution 
as a function of /Cj. For A-^ < 0.1 wehavee max = 0.5 (correspond- 
ing to the plateau in Fig. |4j». For higher values of acJ this function 
is well fitted by a Pade approximant of the form 



2 m =0 ^m^s)'" 



(B.l) 



where a = 0.502, a x = -0.301, a 2 = 0.043, a 3 = 0.078, 
a A = -0.037 and b = 0.932, b\ = 0.092, b 2 = -0.107, which 
provides an excellent fit (x 2 < 4 x 10~ 6 ), for values of in the 
range [0.1,1.5]. 

B.2. Fitting formulae for the ellipticity of the mass distribution 
of the PNFW model 

We have verified that the ratio ss/e is well-fitted by a third-order 
polynomial in s for the whole range of this parameter, 



(B.2) 



The values of the coefficients c, are obtained from this fit for each 
Kj in the considered range. These functions are in turn fitted by 
Pade approximants of the form 



_ 2 n= dn(K^)" 



iJm=0 c m\ K s ) 



(B.3) 



where the coefficients d„ and e m are given in Table IB. ll 





Co(Kj) 




C2(^) 


C3<>j) 


do 


2.523 


-0.380 


0.701 


-0.005 


di 


0.978 


-0.191 


-1.384 


-0.014 


d 2 


2.503 


-2.439 


3.333 


0.110 


d 3 


0.303 


-0.540 


0.188 


-0.217 


d 4 


0.778 




-0.699 




eo 


1.281 


0.183 


0.585 


0.033 


e\ 


0.492 


0.329 


-0.360 


-0.079 


ei 


0.208 


0.669 


1.684 


0.361 


e 3 


0.859 




-0.773 




x 1 


7.9 x 10- 6 


2.8 x 10-" 


1.4 x 


6.7 x 10-' 



Table B.l. Results from the regression analysis using the Pade 
approximant for c,-(/Cj ). Polynomials of different degrees are used 
for co/c 2 and ci/c^. The last row corresponds to the values of^- 2 
for each function c,(Kj). 



We verified that the combination of lB.2l and IB.3l is indeed a 
good approximation for the ellipticity of the convergence con- 
tours in the region associated to the curves = +R t h- For ex- 
ample, the absolute difference £e(s,/Cj) - e^ f (e, k^R^ - +10)| 
is at most 10 3 in the whole range of s and k v s . 

B.3. Fitting formulae for mapping characteristic 
convergences 

The relation between the characteristic convergence of the 
PNFW model (k 1 ^) and the characteristic convergence of the 
ENFW model («f ) is fitted by a quadratic function in k v s 



i^is,^) = p (s) + pds)^ + P2(e)(^) 2 



(B.4) 



By construction po(s) — > 0, pi(s) — > 1, and p 2 (_s) — » as s — > 0. 
We divide the regression analysis into two ranges of . For 
< 0.1 we choose po{s) = and found that p\(s) and p 2 (s) are 

well-fitted by 



Pi(e) = ^8nS n , p 2 (e) 



n=0 



ho + h\E + h 2 e 
ko + k\8 



(B.5) 



with coefficients g„, h„ and k„ given in Table IB.2I (which are 
valid for s < 0.8). We check the accuracy of Eq. (IB.4t with 
p\(e) and p 2 (s) given above by computing the x 2 for a* < 0.1 
and s < 0.6, and we found that it is less than 3 x 10~ 9 in this 
range of parameter values. 





Pi( £ ) 






go 


1.00 


ho 


0.001 


Si 


0.005 


h 


0.446 


82 


-0.035 




-0.039 


83 


0.316 




0.905 


g4 


-0.257 


h 


-0.845 


2 

X' 


1.25 x 10- 8 




2.56 x 10-' 



Table B.2. Results from the regression analysis using a polyno- 
mial form for p\(s) and a Pade approximant for p 2 (s). The last 
row corresponds to the values of^ 2 for each function. 



For 0.1 < *Cj < 1.5 the functions po(s) and p 2 (e) are well 
fitted by polynomials 



Po,2(e) = ^ q„e n , 



(B.6) 



with coefficients q n given in Table lB.3l whereas p\(s) is fitted by 

1 + S\S + S2S 2 + S3£ 3 



1 + t\S + t 2 E 2 



(B.7) 



with si = -0.353,^2 = -0.0270, s 3 = -0.133 and t\ = 
-0.339, t 2 = -0.473, which gives x 2 = 8.55 x 10~ 7 for values of 
e < 0.8 in the considered range. 

We have also checked the accuracy of Eq. dB.4b with 
Po(s), p\(e), and p 2 (s) given above for values of the 
convergence close to the critical curves. We computed 
|«f (e, *f ) - a-j (e, /Cj , R A - ±10)| and found that it is at most 10~ 2 
for e < 0.6. The x 2 of the fit is less than 9 x 10~ 5 in the same 
parameter range. 
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Po( E ) 


P2( E ) 


q\ 


0.008 


-0.030 


12 


-0.091 


0.596 


93 


0.102 


-0.970 


94 


-0.052 


0.787 


95 


-0.034 


0.056 


x 1 


1.25 x 10-' 


1. x 10-" 



Table B.3. Results from the regression analysis using polyno- 
mial forms for po and p2- The last row corresponds to the values 
of x 2 for each function. 
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